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Right-angled hexagon tilings of the hyperbolic plane 


R. Kenyon* 


1 Introduction 

In the spectrum of random planar structures studied in mathematics, the most basic exam¬ 
ples are lattice-based systems like random tilings of the Euclidean plane with polyominoes; 
in these the randomness is encoded in the combinatorial arrangement of tiles. A similar 
basic setting is in classical lattice statistical mechanics, where the randomness is encoded in 
the states of “particles” sitting at lattice sites (typically with a hnite number of states), and 
interacting with nearest neighbors. 

At the other end of the spectrum of random planar processes, one can consider ran¬ 
dom metrics on as in the case of Liouville quantum gravity, where there are no local 
interactions or combinatorics. 

Between these two extremes of either pure combinatorics and fixed geometry, or no 
combinatorics and random geometry, is the case of non-lattice statistical mechanics, in which 
particles are free to change position but still interact with each other. Very few examples of 
this type have been successfully studied. 

We discuss here a case of this intermediate type, where the randomness is encoded in the 
spatial positions of interacting components. 

We work in the hyperbolic plane H, and our components are right-angled hexagons 
(RAHs) of varying shape. The components “interact” in such a way as to form global, edge- 
to-edge tilings of the hyperbolic plane. One can take a different point of view on the same 
system and consider the individual components to be bi-infinite geodesics, which interact 
with each other so as to intersect orthogonally (and so that the complementary components 
are hexagons). 

Other models of random geometric structures in have been studied, see e.g. [2]. 
However in these cases the underlying measures are simpler in a Markovian sense: once a 
geodesic is determined the left and right half-space structures are independent. Although 
our model eventually boils down to a similar argument, the analogous condition is a priori 
much less evident. 
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There is a periodic tiling with regular (all sides of equal length) RAHs, see Figure]^ It 
is not at all clear that this structure is flexible in a way that preserve angles. 



Figure 1: The tiling with regular RAHs. 

Let 0 G El denote the center of the Poincare disk. Let be the space of edge-to-edge 
tilings of the hyperbolic plane El with RAHs. It has a natural topology: two tilings Ti, T 2 are 
close if there is a homeomorphism 0 : El —)■ El which is £-close to the identity {\4>{z) — z\ < e) 
on a large neighborhood of 0. The space also comes with a natural action of the isometry 
group of El; it is thus a Riemann surface lamination, whose leaves are generically hyperbolic 
planes. 

One way to explore O is to study invariant measures on it. Given a measure /i on RAHs 
one can try to construct a measure on O for which the marginal (the induced measure on 
tiles) is given by /i. We prove this for a family of measures p: 

Theorem 1. Given a probability measure v on M_|_ of full support, let /i be the measure on 
RAHs whose edge lengths £5 are independent and distributed according to v. Then there 
is an IsomfM)-invariant probability measure on O of full support in which the tile containing 
the origin (and therefore every tile) has measure pi. 

An analogous result holds for tilings with right-angled pentagons (Section and pre¬ 
sumably for fc-gons with k > 6 as well. For shapes with fewer sides, such as quadrilaterals 
with angles ‘Inj'b, say, we have been unable to formulate a similar result, probably because 
the dimension of the space of shapes is too small. 

One way to state the invariance under Isom(EI) of a measure on O is as follows. Take a 
p-random tiling. Pick an isometry 0 of El independent of the tiling. Applying cf to the tiling 
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results in a new /i-random tiling. In fact, since all tiles have the same area, the random tiling 
“looks the same around any tile”, in the following sense. Choose from the tiling any tile in a 
combinatorial manner, that is, which does not depend on the geometry of the current tiling, 
except for the information about which tile contains the origin. Translate this tile so that 
the origin 0 G El is at a uniform random location in it, and then choose a random rotation 
hxing the origin. Then the result is another exact sample of the measure. 

There are leaves of hi which are closed Riemann surfaces, corresponding to RAH tilings 
which are periodic under the action of a cocompact subgroup of isometries. These leaves come 
in families whose dimension is (as one can show) half the dimension of the Teichmuller space 
of the associated Riemann surface. Each such a family supports many isometry-invariant 
measures, supported on group-invariant RAH-tilings. We will not address here the question 
of what natural measures are supported on these subspaces; we are concerned rather with 
measures of full support in H. 

Another motivation for studying the space H comes from integrable systems. There is a 
surprising connection between the geometry of an RAH and the Yang-Baxter equation for 
the Ising model. The Yang-Baxter equation, or star-triangle equation, for the Ising model 
is a local rearrangement of a graph preserving the so-called partition function of the model. 
The Yang-Baxter equation is a hallmark of an underlying integrable system, and for the Ising 
model one which has been studied in [1]. One can translate the Ising integrable structure 
into an integrable structure on certain spaces of surfaces tiled by RAHs. The measures we 
discuss below are relevant to this system: the “trigonometric” configurations here are in fact 
hxed points of the system, and the measure fiQ discussed below is an invariant measure. This 
will be the subject of a forthcoming work. 

Acknowledgements. We thank Oded Schramm and Andrei Okounkov for helpful discus¬ 
sions, and the referee for several suggestions for improvments. 

2 Measures on RAHs 

2.1 Hexagons 

The side lengths (^i)i=i,,..,6 of a right-angled hexagon satisfy the following relation. Let 
(a, R, c. A, b, C) = (e^b ■ ■ • 5 If c are known (but arbitrary in (0, cxd)) then the re- 
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maining three are determined by 


A + 1 
A - 1 


B + l 
B - 1 


C + l 
C -I 


(1 + abc){a + be) 
{b + ac)(c + ab) 

^(1 + abc){b + ac) 
(a + bc){c + ab) 

^(1 + abc){c + ab) 
(a + bc){b + ac) 


( 1 ) 


It is straightforward to see that the map 'h: (a, b, c) h->• [A, B, C) is an involution. In partic¬ 
ular the family of RAHs is homeomorphic to parametrized by three non-adjacent edge 
lengths. 

One particularly nice family of hexagons are the trigonometric ones, where 


1 — a — b — c — ab — ac — ab + abc = 0, 


or equivalently 

a = tana, 6 = tan/ 3 , c = tan7, 

and a-|-/3-l-7 = |. In this case opposite sides have equal lengths {A = a, B = b, C = c, as 
one can check from ([^). One can dehne a measure on trigonometric RAHs by choosing, for 
example, a, ( 3 ,7 uniformly with respect to Lebesgue measure on the simplex {a-f/3-f7 = |}. 

A more interesting measure /io, of full support in K+, is given by choosing £1,^3, £5 
independently and each distributed with respect to the measure z/q on (0, 00) with density 


F{i) di 


Cdi 

sinh^/3(£) 


where 

P V3r(|)r(|) 

27r3/2 

In terms of u = G ( 1 , cxd) the density is 


f{u) du 


2^/^Cdu 

„ 2 / 3(„2 _ 1 ) 1 / 3 - 


( 2 ) 


( 3 ) 


The truly remarkable property of this measure is 


Theorem 2. If are chosen i.i.d. with distribution vq then £2, 0,^6 will be i.i.d. with 

distribution uq as well. 
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Proof. This is a computation: one needs to show that 


/(a)/(6)/(c) da db dc = f{A)f{B)f{C) dA dB dC 


that is, that the Jacobian of the mapping from (a, b, c) to (A, B, C) defined by Q is 



/(q)/W/(c) 

f{A)f{B)f{Cy 


□ 


One can “discover” the distribution po as follows. If we assume a, 6, c are i.i.d. with some 
density g{x)dx then A,B,C will be distributed according to the density 



^*{g{a)g{b)g{c)dadbdc) = g{a)g{b)g{c) 


Massaging the right-hand side of this expression into a form separating out the A, B and 
C dependence, one sees that there is a unique choice of g (up to scale), given by ([^, for 
which it can be written as g{A)g{B)g{C)dAdB dC. It follows that po is fhe only probability 
measure with this “independence-preserving” property. 

Question. Does this probability measure /io have a geometric signihcance? 

One can of course choose any other measure u on IR+ and dehne a corresponding measure 
/X = niy) on RAHs by choosing lengths £1, £3, £5 independent and ix-distributed. The remain¬ 
ing lengths are then determined (and will have the same marginal distribution as 

each other but will not be independent and typically not distributed according to v). The 
resulting measure /x = fi{u) on RAHs is invariant under rotation by two indices. 


3 Tiling 


3.1 Hextrees 

We will build an Isom(]HI)-invariant probability measure A on D with the property that the 
marginal distribution of a tile (the probability measure restricted to any single tile) is /x(z^) 
for some u as above. The building blocks of this measure are collections of RAHs glued 
together in a ternary-tree structure as in Figure We call such sets hextrees. That is, 
for each edge of the inhnite regular degree -3 tree take an independent random real variable 
distributed according to v. Then, for each vertex of the tree, use the three adjacent random 
variables as lengths of edges 1 , 3,5 of a RAH. These RAHs are glued together so that RAHs 
at adjacent vertices are glued along the corresponding edge. We refer to these glued edges 
as the cross-edges of the hextree, since they cross between boundary components. 
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Figure 2: A random hextree in H. For the purposes of illustration this one was constructed 
using a measure on side lengths which is more tame than i/q- 


3.2 Gluing 

Now the construction of the tiling in hi is as follows. Start with a hextree as above; for each 
of its boundary geodesics construct another hextree, conditioned on having a boundary with 
the same subdivision points as the hrst, but otherwise independent of the hrst (we discuss 
how this is accomplished below); then glue these hextrees along their common boundary 
so that the individual RAHs meet edge-to-edge. We continue gluing hextrees to existing 
boundaries in this way, always choosing independent new hextrees conditioned to agree 
along the common boundary. 

To show that the result is a tiling in hi, it remains to show all of El is covered. This is 


proved in section |3.3| below. 

How does one construct a hextree with the same boundary data (along one boundary 
geodesic) as a given hextree? Given a hextree R let Y be the boundary data (the sequence 
of edge lengths) along one boundary of i?, and X be the remaining lengths dehning the rest 
of the hextree. 

There is a general theorem which applies in this situation [3]: given a joint random 
variable {X,Y) G Gi x (22 (that is, a random variable in a product space), there is a 
conditional probability measure {X,Y)\Y obtained by conditioning on T; a random sample 
(A', Y)\Y of this measure can be chosen independently of (X, Y) conditional on Y. So to glue 
a hextree to another we just resample this conditional measure. In practice this resampling 


is complicated; see section |3.4| below. 

The Mobius invariance of the resulting measure is a consequence of the reversibility of 
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Figure 3: A RAH tiling. For this figure we used the measure v on edge lengths with density 
(for u = e^) f{u) = C{u — with mean 4. 

the gluing procedure: each new hextree is distributed according to the same measure /i, 
and the joint probability measure on the glued object is symmetric: one cannot tell where 
the construction began and in what order the hextrees were glued. Furthermore, within 
each hextree each RAH has the same distribution, with edges 1, 3, 5 which are independent, 
so there is no information about which RAH was the “first”. In particular, as mentioned 
earlier, if we choose any RAH in the tiling translate it so that the origin 0 G El is at a uniform 
random location in that RAH, the result is another exact sample of A. 

Question. Is the measure on tilings constructed from the above measure yUo characterized 
by some natural property like entropy maximization? 

Question. Is there an Isom(]HI)-invariant measure on H supported on trigonometric RAHs? 

3.3 Covering property 

To show that the union of hextrees covers all of El, we need to show that hextrees are typically 
“thick”: if some cross-edge of a hextree is short then along the geodesic containing this short 
cross-edge the other cross-edges of nearby hextrees are typically no shorter. 

This follows from the reversibility of the construction of the tiling. Let xq be a cross-edge. 
Choose one of the two cross-edges adjacent to it on the same geodesic. This cross-edge xi 
is just as likely to be greater than xq as smaller than xq, by reversibility. Similarly for the 
next cross-edge X 2 , and so on. So the chance of getting a sequence of smaller and smaller 
cross-edges starting from xq goes to zero with the length of this sequence. 
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Thus along the “cross-edge geodesics” there are no accumulation points of vertices of 
RAHs, and so all of El is covered. 

3.4 Resampling 

We are given a hextree R with boundary data Y = ?/_i, j/q, yi,... }, and wish to con¬ 

struct another, independent of R conditional on having the same boundary data. It suffices 
to construct only the layer of RAHs adjacent to the boundary, since the remaining RAHs can 
be chosen by choosing their cross-edge lengths independently of everything as in the original 
definition of p. To construct the RAHs adjacent to the boundary H, it suffices to construct 
the sequence Xi of edge lengths of the edges incident to Y (ending on Y and perpendicular 
to it, so that Xi ends between yt and yi+i) since the three consecutive edges Xi-i,yi,Xi will 
determine the RAH at position i. 

We construct the Xi using a successive approximation procedure, starting initially from 
the corresponding lengths x[ of R, that is, from the edge lengths of the edges of R incident 
to Y from the other side. 

We use the fact that the sequence {(xj,?/*)} is Markovian: given one pair {xi,yi) the 
values {xj,yj) for j > i are independent of the values {xj,yj) for j < i. This follows from the 
construction of fi from the sequence of independent cross-edge lengths: the values {xj, yj) for 
j > i depend only on the choices of cross-edge lengths on one side of Xj, the values {xj,yj) 
for j < i depend only on the cross-edge lengths on the other side of x*. 

From the sequence {x'jjgz, erase all X 2 i for all i. Then for each z, resample X 2 j according 
to the conditional measure defined by x' 2 i_i, ?/ 2 i, a: 2 i+i- By this we mean, by the Markov 

property, the marginal distribution of x' 2 i (given everything else) only depends on these four 
values. Resampling each X 2 j is then a straightforward finite-dimensional computation. 

Now repeat for the odd indices: erase each X 2 j_,_i and resample it according to its four 
neighboring values. Now iterate, resampling successively the even values then the odd values 
over and over. The resampling map preserves the measure y conditioned on Y, and the Xj’s 
will quickly become decorrelated from their initial values, that is, after a small number of 
iterations the state is independent of the initial state. To be a little more precise, if we are 
interested in functions of {xj} which depend only on indices in the interval [—m,m], then 
after a polynomial number (in m) of iterations the system is close to uncorrelated. These 
facts have been established in [1] in a similar setting. 

Although it should be possible to show that the dependence of x, on yj decays exponen¬ 
tially with \i — j\, we have not attempted to prove this here. 


4 Pentagons 

Let £i,... ,£5 be the side lengths of a right-angled pentagon (RAP). Any two of these deter¬ 
mine the rest: we have for example 


cosh £4 = sinh^i sinhf'2, 

and its cyclic rotations. Explicitly, given £1,^25 which must satisfy sinh£isinh£2 > 1 , we 
have 

cosh £1 

a/ sinh^ tx sinh^ £2 — 1 

a/ sinh^ £1 sinh^ £2 — 1 
cosh ^2 

a/ sinh^ £1 sinh^ £2 — 1 


sinh £3 = 

sinh £4 = 
sinh £5 = 


There is a natural measure r on right-angled pentagons (RAPs) which is invariant under 
rotations, dehned as follows (one can discover this measure in a similar manner as one 
discovers the measure vq as discussed after the proof of Theorem |^ . If we parametrize by 
£1, £2, fhe density is 


/(^ 1 ; -^2) df'i M2 


CM1M2 

\/ sinh^ ii sinh^ £2 — 1 


CM1M2 

sinh £4 


where C = iP(l/A/ 2 )^, K being the complete elliptic integral of the hrst kind. For this 
measure each edge has expected length n/2. 

In terms of exponential of the lengths we have (with a = and 6 = 


f{a, b) da db 


da db 

a/ {ab -|- a -|- 6 — 1) {ab — a -|- 6 -|- 1) {ab -|- a — 6 -|- 1) {ab — a — b — 1) 


As in the case of right-angled hexagons, one can construct an Isom(]HI)-invariant measure 
on tilings of El with RAPs whose marginal distribution is r, as follows. 

Gluing four RAPs edge-to-edge around a vertex produces a right-angled octagon after 
merging four pairs of edges. The lengths of its four edge lengths opposite to the 

central vertex (which were originally edges of pentagons) satisfy one relation: 


cosh cosh £5 = cosh £3 cosh G- 


( 4 ) 


Given four lengths satisfying this relation, there is a one-parameter family of choices for the 
central four pentagon edges. 

We can glue such octagons together along the G for i odd, choosing lengths at each 
glued edge which are independent subject to (|^— as explained below—to make a thickened 
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4-valent tree, an “octatree”. Then we can proceed as above, gluing octatrees to make a tiling 
of e. 

To choose the edge lengths of the cross-edges of the 4-valent tree, subject to Q at each 
vertex, proceed as follows. Take a bi-in£nite straight path (not turning left or right at a 
vertex) in the tree through a hxed vertex n; assign values to its edges which are i.i.d. and 
distributed according to the cosh of the edge lengths. Assign the edges along the other 
straight path through v also i.i.d. except that the values on the two edges adjacent to v must 
have product equal to the product of values on the other two edges adjacent to v. Now 
proceed outwardly from n; along each new bi-in£nite geodesic encountered, assign values on 
it i.i.d. subject to the condition (|^ at each vertex. Since this is a tree the reversibility of the 
assignments guarantees that the resulting object is Isom(]HI)-invariant. 
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